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Circularly-polarized  extreme  ultraviolet  and  X-ray  radiation  is  useful  for  analysing  the  structural,  electronic  and  magnetic 
properties  of  materials.  To  date,  such  radiation  has  only  been  available  at  large-scale  X-ray  facilities  such  as  synchrotrons. 
Here,  we  demonstrate  the  first  bright,  phase-matched,  extreme  ultraviolet  circularly-polarized  high  harmonics  source.  The 
harmonics  are  emitted  when  bi-chromatic  counter-rotating  circularly-polarized  laser  pulses  field-ionize  a  gas  in  a  hollow- 
core  waveguide.  We  use  this  new  light  source  for  magnetic  circular  dichroism  measurements  at  the  /W-shell  absorption 
edges  of  Co.  We  show  that  phase-matching  of  circularly-polarized  harmonics  is  unique  and  robust,  producing  a  photon  flux 
comparable  to  linearly  polarized  high  harmonic  sources.  This  work  represents  a  critical  advance  towards  the  development 
of  table-top  systems  for  element-specific  imaging  and  spectroscopy  of  multiple  elements  simultaneously  in  magnetic  and 
other  chiral  media  with  very  high  spatial  and  temporal  resolution. 


Circularly-polarized  radiation  in  the  extreme  ultraviolet  (EUV) 
and  soft  X-ray  spectral  regions  has  proven  to  be  extremely 
useful  for  investigating  chirality- sensitive  light-matter  inter¬ 
actions.  It  enables  studies  of  chiral  molecules  using  photoelectron 
circular  dichroism1,  ultrafast  molecular  decay  dynamics2,  the 
direct  measurement  of  quantum  phases  (for  example,  Berry’s 
phase  and  pseudo -spin)  in  graphene  and  topological  insulators3,4 
and  reconstruction  of  band  structure  and  modal  phases  in  solids5. 
For  magnetic  materials,  circularly-polarized  soft  X-rays  are  particu¬ 
larly  useful  for  X-ray  magnetic  circular  dichroism  (XMCD)  spec¬ 
troscopy6,  which  enables  element-selective  probing  as  well  as 
coherent  imaging  and  holography  of  magnetic  structures  with  nano¬ 
metre  resolution7-9.  XMCD  can  also  be  used  to  extract  detailed 
information  about  the  magnetic  state  by  distinguishing  between 
the  spin  and  orbital  magnetic  moments  of  each  element.  Thus, 
time- resolved  XMCD  can  probe  the  element- specific  dynamics  of 
the  spin  and  orbital  moments  when  interacting  with  the  electronic 
and  phononic  degrees  of  freedom  in  a  material10,11.  Although  ultra¬ 
fast  magnetic  dynamics  has  been  studied  for  20  years12,  the  time 
resolution  of  XMCD  remains  >100  fs,  limited  by  the  pulse  duration 
and  timing  jitter  of  synchrotron  and  free-electron  laser  pulses13,14. 
To  date,  it  has  not  been  possible  to  simultaneously  probe  the 
spin  dynamics  of  multiple  elements  within  a  single  sample, 
because  the  photon  energy  must  be  tuned  across  the  various  absorp¬ 
tion  edges  at  large-scale  facilities  (where  these  experiments  are 
currently  performed). 

Table-top  soft  X-ray  sources  based  on  high  harmonic  upconver- 
sion  of  femtosecond  laser  pulses  represent  a  viable  alternative  to 
large-scale  sources  for  many  applications,  due  to  their  unique 
ability  to  generate  bright,  broadband,  ultrashort  and  coherent 
light  with  an  energy  spectrum  reaching  into  the  keV  region15. 
High  harmonic  generation  (HHG)  not  only  enables  coherent 
imaging  of  nanometre  structures  with  a  spatial  resolution  approach¬ 
ing  the  diffraction  limit16,  but  also  accesses  the  fastest  dynamics  in 


atoms,  molecules,  solids  and  plasmas  with  unprecedented  time 
resolution17-23.  The  large  bandwidth  and  short  temporal  duration 
of  HHG  makes  it  a  versatile,  element-specific  probe  of  coupled 
spin,  charge  and  phonon  dynamics  on  ultrafast  timescales.  In  the 
case  of  magnetic  materials,  recent  work  has  shown  that  HHG  can 
simultaneously  probe  the  magnetic  state  of  the  3d  ferromagnets 
Fe,  Co  and  Ni  by  taking  advantage  of  the  fact  that  the  reflectivity 
near  the  M-shell  absorption  edges  depends  on  the  orientation  as 
well  as  the  magnitude  of  the  magnetization24.  This  novel  capability 
has  made  it  possible  to  gain  a  new  fundamental  understanding 
about  the  timescale  of  the  exchange  interaction,  as  well  as  revealing 
the  important  roles  of  local  spin  scattering  and  non-local  spin  trans¬ 
port,  all  of  which  can  occur  on  few-femtosecond  timescales25,26.  The 
vast  majority  of  HHG  applications  to  date,  however,  have  used  lin¬ 
early  polarized  light,  which  can  be  generated  relatively  efficiently. 

For  decades  it  was  generally  assumed  that  HHG  from  atoms  is 
bright  when  both  the  driving  laser  and  HHG  fields  are  linearly  polar¬ 
ized.  The  source  of  this  misconception  is  that  HHG  is  a  re-collision 
phenomenon.  A  bound  electron  is  ripped  from  an  atom  by  the 
intense  electric  field  of  a  driving  laser  pulse  and  then  accelerated  as 
a  free  electron  until  it  re-encounters  its  parent  ion27.  For  a  linearly 
polarized  driving  field,  the  electron  accelerates  on  a  linear  trajectory, 
and  therefore  easily  re-collides  with  the  parent  ion.  When  driven  by 
a  slightly  elliptical  laser  field,  the  electron  has  some  probability  of 
re-colliding  with  its  parent  ion  due  to  lateral  spreading  (quantum 
diffusion)  of  the  wavefunction  in  the  continuum28.  This  results  in 
the  generation  of  slightly  elliptically  polarized  high  harmonics29. 
In  contrast,  for  circularly-polarized  driving  lasers  (or  elliptically- 
polarized  lasers  with  a  large  ellipticity),  the  probability  of  re- collision 
and  the  emission  of  high  harmonics  is  completely  suppressed. 
Nevertheless,  because  of  important  applications  in  materials  science, 
there  was  strong  motivation  to  generate  bright  circularly-polarized 
high  harmonic  beams30-40.  The  first  experimentally  measured 
circularly-polarized  high-order  harmonics  were  produced  by 
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converting  the  polarization  of  linearly  polarized  HHG  to  circular, 
using  a  reflective  quarter-waveplate33.  This  approach  is  very  lossy,  is 
limited  to  narrow  spectral  regions  constrained  by  the  available  multi¬ 
layer  mirror  materials,  and  to  date  has  not  proven  useful  for  appli¬ 
cations.  A  direct  approach  for  generating  circularly-polarized  HHG 
was  suggested  almost  two  decades  ago41,42  and  recently  measured  by 
Fleischer  and  colleagues43.  In  this  scheme,  circularly-polarized  high 
harmonics  are  driven  by  co-propagating  circularly-polarized  bi-chro- 
matic  fields  that  rotate  in  opposite  directions  and  interact  with  an  iso¬ 
tropic  gas.  Notably,  Fleischer  and  co-workers  measured  the  ellipticity 
of  weak  harmonics  generated  from  HHG  in  a  thin  gas  jet,  but  not  their 
helicity  (that  is,  right-  versus  left-circular  polarization). 

In  this  Article  we  generate  bright,  phase-matched,  circularly- 
polarized  high  harmonic  beams  for  the  first  time,  and  then  use  this 
unique  light  source  to  implement  the  first  table-top  XMCD  measure¬ 
ments.  To  achieve  this,  we  first  demonstrate  (both  experimentally 
and  theoretically)  that  phase-matching  of  high  harmonics  driven 
by  circularly  polarized  bi- chromatic  waves  differs  significantly 
from  all  other  approaches  to  phase-matching  to  date.  Phase¬ 
matching  is  more  robust  for  circularly-polarized  bi- chromatic 
drivers  than  for  HHG  with  bi-chromatic  linearly-  or  elliptically- 
polarized  drivers.  Moreover,  phase-matching  can  selectively 
enhance  either  left  or  right  circularly-polarized  harmonics,  a  critical 
feature  for  generating  circularly-polarized  attosecond  pulses. 
Using  this  new  understanding,  we  generate  bright  circularly- 
polarized  high  harmonics  by  co-propagating  bi-chromatic  driving 
laser  beams  (the  fundamental  and  second-harmonic  beams  from  a 
Ti:sapphire  laser)  that  are  circularly-polarized  with  opposite 
helicity  in  a  gas-filled  hollow  waveguide  (see  Supplementary 
Section  6  for  details  on  the  important  features  of  this  source).  We 
show  that  bright,  circularly-polarized,  high  harmonics  in  Ne  can 
span  multiple  M-edges  of  the  3 d  ferromagnets,  and  obtain  a 
photon  flux  of  x  109  HHG  photons  per  second  per  harmonic 
(Supplementary  Section  5),  which  is  comparable  to  the  linearly- 
polarized  HHG  flux  that  was  used  very  successfully  in  previous  ultra¬ 
fast  element- selective  magneto-optic  HHG-based  experiments24-26. 
Finally,  we  measure  the  circular  dichroism  of  Co  throughout  the 
M-shell  absorption  edge  spectral  region.  This  measurement  is  also 
the  first  to  measure  the  helicity  of  circularly-polarized  high-order 
harmonic  beams.  In  the  future,  the  use  of  mid-infrared  pump 
lasers  should  make  it  possible  to  generate  bright  circularly-polarized 
high  harmonics  at  the  L-shell  absorption  edges  of  Fe,  Co  and  Ni,  in  a 
photon  energy  range  approaching  1  keV  (ref.  15). 

Phase-matching  of  circularly-polarized  HHG 

We  first  describe  the  theoretical  spectral  and  polarization  features  of 
our  source  driven  by  a  left-rotating  circularly-polarized  fundamental 
beam  at  a  central  wavelength  of  790  nm,  and  a  right- rotating  circu¬ 
larly-polarized  second -harmonic  beam  (Fig.  la).  We  consider  an 
idealized  case  in  which  the  medium  is  isotropic  and  time-independent, 
and  where  the  driving  laser  is  perfectly  periodic.  Independently  of  the 
relative  intensity  or  phase  between  the  two  pump  fields,  the  bi-chro- 
matic  field,  E BC(f),  exhibits  the  following  dynamical  symmetry: 
Ebc(^  +  TI3)  =  -R(120°)  EBC(f).  T  is  the  optical  cycle  of  the  fundamental 
laser  and  -R(12o°)  is  the  120°  rotation  operator  in  the  polarization  plane. 
The  emitted  HHG  field  conforms  to  the  same  dynamical  symmetry30. 
Thus,  as  can  be  verified  by  taking  the  Fourier  transform  of  the  dyna¬ 
mical  symmetry  constraint,  the  complex  amplitude  of  each  ^th- order 
harmonic,  Eq,  satisfies  the  following  eigenvalue  equation 

Eqe-2m"=R(mo)Eq  (i) 

The  solutions  of  equation  (1)  are  a  left  circularly-polarized  field 
with  eigenvalue  e~2ml 3  and  a  right  circularly-polarized  field  with 
eigenvalue  e+2ml3.  The  left  circularly-polarized  gth-order  harmonic, 

100 


Figure  1  |  Pump  symmetry  and  phase-matching  condition  for 
circularly-polarized  HHG.  a,  The  combined  electric  field  of  a  left-circularly- 
polarized  790  nm  driving  laser  beam  (red)  and  a  counter-rotating  second- 
harmonic  field  (blue)  is  a  threefold  rosette  shape  (purple).  The  system  has  a 
discrete  rotation  dynamical  symmetry  of  the  light  field,  where  a  delay  of 
7/3  acts  as  a  120°  rotation.  This  dynamical  symmetry  gives  rise  to 
circularly-polarized  harmonics  where  harmonic  orders  q  =  3m  -  1  rotate  right 
and  g  =  3m  +  1  rotate  left,  b,  Calculated  coherence  length  versus  HHG  order 
under  phase-matching  conditions  that  optimize  the  28th  harmonic.  Phase¬ 
matching  of  circularly-polarized  HHG  is  helicity-selective,  so  HHG  orders 
<7  =  3m  + 1  (left-circular  polarization,  LCP)  have  a  longer  coherence  length 
(better  phase-matching)  than  <7  =  3/77-1  (right-circular  polarization,  RCP). 
Reabsorption  of  high  harmonics  in  the  gas  medium  reduces  the  contrast 
arising  from  this  effect,  as  detailed  in  Supplementary  Section  1. 

E^l,  satisfies  e~^2mq/3^EqL  =  e-(2ra/3)E^  L,  which  is  fulfilled  only  when 
q  =  3m+l.  Similarly,  the  right  circularly-polarized  gth-order 
harmonic,  E^R,  satisfies  c-(2tt^/3)E^r  =  e+(2ra/3)E^R,  which  is  fulfilled 
only  when  q-3m-\.  Notably,  q  =  3m  harmonics  do  not  satisfy 
equation  (1)  and  are  therefore  theoretically  forbidden. 

Phase-matching  is  essential  for  obtaining  a  bright  and  collimated 
high  harmonic  beam.  In  HHG  with  linearly-polarized  bi-chromatic 
drivers,  the  phases  and  intensities  of  the  emitted  harmonics  depend 
strongly  on  the  relative  phase  between  the  chromatic  drivers,  due  to 
the  intrinsic  phases  of  the  rescattering  electrons44.  Consequentially, 
full  phase-matching  in  this  case  requires  that  the  two  driver  fields 
and  the  harmonic  field  propagate  at  the  same  phase  velocity,  a 
very  challenging  condition  to  satisfy.  Fortunately,  phase-matching 
of  HHG  driven  by  circularly-polarized  bi-chromatic  drivers  is 
very  different.  In  this  case,  the  relative  phase  between  the  chromatic 
driver  components  does  not  influence  the  shape  (rosette)  of  the  bi- 
chromatic  driver  (see  Supplementary  Movie),  and  hence  it  also  does 
not  influence  the  path  length  of  the  re-colliding  electron  and  the 
intrinsic  phase  it  accumulates.  Thus,  phase  mismatch  of  the 
process  is  determined  only  by  the  extrinsic  phase,  and  is  described 
by  the  ‘ordinary’  phase-mismatch  equation  in  nonlinear  optics: 

—  =  A k  =  kq-lk1-  mk2  (2) 
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Figure  2  |  Experimental  apparatus.  Circularly-polarized  harmonics  are  generated  by  focusing  fundamental  and  second-harmonic  pump  beams  (with 
opposing  helicities)  into  a  gas-filled  hollow  waveguide.  The  polarization  states  of  the  two  drivers  can  be  adjusted  independently,  thereby  controlling  the 
helicity  of  the  generated  harmonics.  The  magnetic  state  of  a  sample  can  then  be  measured  with  elemental  specificity  via  XMCD  by  using  harmonics  that 
overlap  with  an  appropriate  core-level  absorption  edge,  for  example,  the  A/I-edge  of  Co.  Right  inset:  helicity  selectiveness  for  phase-matching  of  circularly- 
polarized  HHG,  where  left-rotating  harmonics  (q  =  3m  +  1)  and  right-rotating  harmonics  (q  =  3m  -1)  cannot  be  phase-matched  simultaneously.  In  other 
words,  if  the  q  =  3m  +  1  harmonics  are  phase-matched  with  the  pump  lasers,  then  the  q  =  3m  -  1  harmonics  are  mismatched  by  Ak.  Left  inset:  typical 
experimental  spectrum  recorded  on  the  CCD. 


where  kq ,  kx  and  k2  are  the  wavevectors  for  the  ^th- order 
high-harmonic,  fundamental  and  second  harmonic  beams,  respect¬ 
ively,  A k  is  the  phase  mismatch  of  the  upconversion  process,  and  Lc 
is  the  resulting  coherence  length,  i  and  m  are  the  number  of 
photons  of  the  fundamental  and  second-harmonic  pump  beams, 
respectively,  converted  to  a  gth-order  high  harmonic  photon. 
Conservation  of  energy  and  angular  momentum  impose  that 
q  =  i  +  2m,  and  i  =  m  ±  1,  respectively43.  These  conservations 
imply  that  left  (right)  circularly-polarized  harmonics  are  generated 
from  a  combination  of  m  second-harmonic  photons  and  m  +  1 
(m  -  1)  fundamental  photons. 

For  simplicity,  we  associate  kx  and  k2  with  the  effective 
index  change,  A  nl  and  An2,  respectively,  defined  as 
A ni  =  (ki/[2n/Xi\  -  1),  where  Af  is  the  vacuum  wavelength  of  the  fth 
driver  (i  =  1  for  the  fundamental  pump  and  i  =  2  for  the  second- 
harmonic  pump).  Thus,  we  can  rewrite  equation  (2)  in  the  form: 

—  =  \m(Anl  +  2An2)  ±  AnJ  (3) 

Lc  'h 

where  the  +  and  -  signs  correspond  to  the  q  =  3  m  +  1  and  q  =  3m  -  1 
groups  of  harmonics,  respectively.  The  factor  ±A n1  leads  to  a  large 
difference  in  the  coherence  lengths  between  the  combs  of  left  and 
right  rotating  circular  harmonics.  Physically,  this  difference  reflects 
the  different  ratios  between  the  fundamental  and  second-harmonic 
photons,  (m  +  1  )/m  and  (m  -  1  )/m,  which  participate  in  the  gener¬ 
ation  of  left  and  right  circularly-polarized  harmonics,  respectively. 

Figure  lb  presents  numerical  results  for  the  phase-matching 
of  circularly-polarized  HHG.  The  coherence  lengths  shown  in 
Fig.  lb  are  calculated  for  A nx  and  A n2  that  optimize  phase-matching 
for  the  28th  harmonic.  As  expected,  the  q  =  3m  +  1  harmonics  are 
brighter  (i.e.  better  phase-matched  to  the  bi-chromatic  pump) 
than  the  harmonics  with  opposite  helicity.  Notably,  to  implement 
phase-matched  HHG  frequency  upconversion,  the  ionization  level 
of  the  gas  must  be  below  a  critical  value,  qcr  (ref.  45).  For  HHG 
with  a  circularly  polarized  bi-chromatic  driver,  qcr  has  a  value 
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between  the  critical  ionization  for  HHG  driven  by  the  fundamental 
or  second  harmonic,  separately  (Supplementary  equation  (S.3)). 
This  also  results  in  a  phase-matching  cutoff  photon  energy  that  is 
somewhere  between  the  cutoff  corresponding  to  the  fundamental 
and  second  harmonic,  separately.  Finally,  we  note  that  our  model 
does  not  include  the  HHG  intensity- dependent  phase,  through 
which  the  temporal  envelope  of  the  pump  pulse  affects  the  instan¬ 
taneous  phase  mismatch.  However,  the  dependence  of  the  harmonic 
phase  on  the  intensity  of  the  bi-chromatic  pump  laser  is  weak38. 

Figure  2  presents  a  schematic  of  the  experimental  apparatus.  A 
Ti:sapphire  oscillator  is  used  in  combination  with  a  single-stage 
regenerative  amplifier,  delivering  sub-45  fs  pulses  at  a  repetition 
rate  of  4  kHz,  which  are  centred  at  a  wavelength  of  790  nm  (red), 
with  an  energy  of  2.5  mj  per  pulse.  After  frequency- doubling  the 
laser  beam  in  a  (3-phase  barium  borate  (BBO)  crystal,  the  bi-chro- 
matic  co -propagating  driving  fields  are  separated  by  a  dichroic 
mirror  (DM)  into  two  different  arms  of  a  Mach-Zehnder  inter¬ 
ferometer.  An  optical  delay  stage  in  the  red  arm  compensates  for 
the  relative  time  delay  between  the  two  colours,  and  the  polarization 
of  each  arm  is  fully  controlled  by  a  pair  of  half-  and  quarter- 
wavelength  retardation  plates.  The  red  (1.6  mj/pulse)  and  blue 
(0.43  mj/pulse)  driving  lasers  are  focused  into  a  150-pm-diameter, 
2-cm-long  gas-filled  hollow  waveguide  using  lenses  with  focal 
lengths  of  50  cm  and  75  cm,  respectively.  The  circularly-polarized 
HHG  beam  that  emerges  from  the  waveguide  then  passes  through 
an  A1  filter  to  block  the  pump  laser  beams.  In  the  first  experiment, 
the  HHG  beam  was  spectrally  dispersed  using  a  spectrometer  com¬ 
posed  of  a  toroidal  mirror,  a  laminar  Au  grating  with  a  groove 
density  of  500  lines/mm  and  a  charge-coupled  device  (CCD) 
camera.  In  the  second  experiment,  the  HHG  beam  passed 
through  a  magnetized  100-nm-thick  Co  foil  tilted  by  45°  with 
respect  to  the  harmonics  beam  before  entering  the  spectrometer. 
The  magnetization  of  the  Co  foil  was  controlled  by  an  external 
electromagnet  that  produced  a  magnetic  field  of  -4.8  x  104  A  m”1 
(60  mT)  parallel  to  the  sample  plane,  exceeding  the  1.6  x  104  A  m-1 
(20  mT)  coercive  field  of  the  foil.  The  alternating  polarity  of  the 
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Figure  3  |  Generation  of  bright  circularly-polarized  harmonics,  a,  HHG  spectrum  from  Ar.  b,  HHG  from  N2  at  70  torr  (black)  and  Ne  at  650  torr  (red).  The 
HHG  spectrum  for  Ne  covers  the  Co  A/I-edge  (black  dashed  curve47).  The  circular  polarization  of  the  HHG  is  manifested  by  suppression  of  the  3m 
harmonics  (15, 18,  ...,  white  arrows  in  c  and  d).  c,  Colour-coded  experimental  (c)  and  simulated  (d)  pressure-dependent  spectrograms  of  HHG  with  Ne. 
Each  notch  mark  on  left  axes  indicates  10  kPa.  Colour  coding  in  c  is  a  logarithm  of  the  photoelectron  counts  per  second  on  the  CCD.  The  experimental  and 
theoretical  spectrograms  exhibit  a  known  feature  of  phase-matching  in  Ne:  higher-order  harmonics  are  phase-matched  through  a  narrower  pressure  range 
than  lower-order  harmonics  (dashed  lines  guide  the  eye),  e,  Measured  (solid  blue  line)  and  normalized  simulated  (dashed  red  line)  spectra  (on  a  log  scale) 
with  650  torr  of  Ne. 


magnet  was  synchronized  to  the  spectrometer  readout,  so  that  we 
were  able  to  acquire  the  transmitted  HHG  spectrum  for  two  (‘up’ 
and  ‘down’)  magnetization  directions.  The  normalized  difference 
of  these  two  spectra  constitutes  the  XMCD  asymmetry. 
Positioning  the  ferromagnetic  sample  at  45°  to  the  HHG  beam 
allowed  observation  of  XMCD  of  a  sample  magnetized  in  plane46. 

We  first  discuss  the  brightness,  spectrum  and  polarization  of  the 
EUV  light  produced  by  this  novel  HHG  source.  Figure  3a  presents 
the  observed  HHG  spectrum  when  the  hollow  waveguide  was  filled 
with  Ar  gas  at  a  pressure  of  100  torr  (13  kPa),  with  the  driving  laser 
beams  blocked  by  A1  filters  of  ~800  nm  thickness.  Figure  3b  displays 
the  observed  HHG  spectra  from  650  torr  (87  kPa)  of  Ne  and  70  torr 
(9  kPa)  of  N2  using  400-nm-thick  A1  filters.  Clearly,  the  q  =  3m 
harmonics  are  almost  completely  suppressed  over  the  entire 
observed  HHG  spectra  for  all  gas  species.  This  feature  indicates 
that  the  polarization  of  the  harmonics  is  nearly  circular43. 

Figure  3c  and  3d,  present,  respectively,  experimental  and  simu¬ 
lated  pressure-dependent  spectrograms  of  HHG  emission  from  a 
Ne-filled  hollow  waveguide  (see  Supplementary  Section  3  for 
details).  Two  features  are  worth  mentioning.  First,  the  3m  HHG 
orders  are  suppressed  throughout  the  pressure  range,  reflecting  the 
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single- atom  dynamical  symmetry  described  in  equation  (1). 
Second,  the  intensity  of  the  plateau  harmonics  dramatically  increases 
at  pressures  above  100  torr  (13  kPa),  indicating  that  they  are  phase- 
matched  throughout  this  pressure  range.  However,  the  cutoff  harmo¬ 
nics  are  phase-matched  within  a  narrower  pressure  range,  as  marked 
by  yellow  dashed  lines.  This  effect,  which  also  exists  in  HHG  from 
linearly-polarized  driving  fields45,  results  from  two  complementary 
effects.  First,  a  larger  effective  index  change,  |  A nx|,  |An2|,  is  required 
for  phase-matching  higher-order  HHG.  Thus,  deviation  from  perfect 
phase-matching  in  equation  (2),  say  1%  of  A nh2,  affects  the  absolute 
value  of  the  coherence  length  Lc  more  drastically.  Second,  the  longer 
reabsorption  length  of  higher  HHG  orders  in  Ne47  makes  the  phase- 
matched  throughput  even  more  sensitive  to  imperfect  phase-match¬ 
ing  conditions.  Notably,  reabsorption  of  the  high  harmonics  in  Ne 
limits  the  distinction  between  full  and  partial  phase- matching. 
Therefore,  the  output  flux  contrast  between  one  perfectly  phase- 
matched  harmonic  (e.g.  28th  order  in  Fig.  lb)  and  the  rest  of  the 
3  m  +  1  harmonic  orders  is  insignificant.  Most  importantly,  the 
current  flux  of  circularly-polarized  photons  produced  by  our 
scheme  is  sufficient  for  application  experiments  such  as  XMCD,  as 
described  in  the  following  section. 
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Figure  4  |  Co  XMCD.  a,  Spectra  of  circularly-polarized  HHG  transmitted  through  100  nm  of  magnetized  Co  with  the  magnetization  direction  either  'up'  (/up, 
blue)  or  'down'  (/down,  red).  The  dashed  black  envelope  shows  the  transmission  spectrum  of  Co  (units  along  the  right  axis)47.  Inset:  difference  in  signal 
between  'up'  and  'down'  magnetization  for  the  34th  and  35th  orders,  including  experimental  error.  Round  arrows  indicate  HHG  helicity.  b,  Experimental 
XMCD  asymmetry  defined  by  (/up  -  /down)/(/uP  +  /down)  alternating  asymmetry  for  adjacent  harmonics  is  consistent  with  the  prediction  of  alternating 
harmonic  helicity.  c,  Magneto-optical  (MO)  absorption  coefficient  A/?,  derived  from  XMCD  asymmetry  of  the  3m  + 1  (red  diamonds)  and  3m  -  1  (blue 
squares)  harmonics.  The  experimental  MO  coefficient  derived  from  the  3m  +  1  harmonics  accurately  retrieves  literature  MO  values  (dashed  line)  without 
fitting  any  parameters.  The  reduced  intensities  and  XMCD  asymmetry  of  the  3m  - 1  harmonics  suggests  that  these  harmonics  have  lower  ellipticity. 


Magnetic  circular  dichroism 

To  demonstrate  the  utility  of  this  table-top  HHG  source  with  circu¬ 
larly-polarized  EUV  light  for  probing  magnetization,  we  measured 
XMCD  of  a  freestanding  thin  Co  foil.  The  imaginary  (that  is, 
absorptive)  part  of  the  refractive  index  n  for  light  propagating  in  a 
magnetically  saturated  film  is  given  by  Im [n]  =  /3  ±  A/3  where  /3  is 
the  average  absorption  coefficient  and  A/3  is  the  dichroic  absorption 
coefficient  for  circularly-polarized  light,  where  the  (±)  marks  the 
absorption  difference  for  left-  versus  right- circularly  polarized 
EUV  photons46,48.  The  dichroic  absorption  scales  linearly  with  the 
component  of  the  magnetization  parallel  to  the  direction  of  light 
propagation.  For  example,  the  absorption  coefficient  of  a  sample 
with  magnetization  M  that  is  below  its  saturation  value  Msat,  point¬ 
ing  at  angle  6  with  respect  to  the  beam  direction,  is 
(M/Msat)  cos  6  •  A/3.  However,  XMCD  requires  a  bright  circularly- 
polarized  soft  X-ray  source,  because  the  non- dichroic  absorption 
coefficient  is  approximately  ten  times  larger  than  the  dichroic  coef¬ 
ficient,  leading  to  strong  average  absorption  that  far  exceeds  the 
dichroic  signal. 

We  measured  the  HHG  spectrum  after  passing  through  a  mag¬ 
netized  Co  foil,  as  shown  in  Figs  2  and  4.  Figure  4a  presents  the 
transmitted  spectrum,  Iup  and  Idown,  for  ‘up’  and  ‘down’  Co  magne¬ 
tization,  respectively.  The  normalized  XMCD  asymmetry, 
A  =  (Iup  -  Idown)  /  (Iup  +  Idown)  is  shown  in  Fig.  4b.  The  asymmetry 
of  the  3m  +  1  and  3m  -  1  harmonics  exhibits  opposite  signs,  proving 
that  the  helicity  of  the  3  m  +  1  harmonics  is  indeed  opposite  to  that 
of  the  3m  -  1  harmonics. 

Finally,  we  used  the  measured  XMCD  asymmetry  to  extract  the 
magneto -optical  (MO)  dichroic  absorption  coefficient,  A/3,  for  Co. 
The  MO  dichroic  absorption  coefficient  is  given  by 
A/3  =  tanh-1(A)/(2|k0|  cos  6L),  where  k0  is  the  light  wavevector,  L 
is  the  optical  path  length  in  the  sample  and  6  is  the  angle  between 
the  sample  magnetization  and  beam  propagation  direction46. 
Because  the  sample  magnetization  is  tilted  by  45°,  the  optical  path 
is  \fl  times  the  sample  thickness  d,  and  cos  6  =  l/\/2>  the  MO 
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dichroic  coefficient  is  simply  A/3  =  tanh_1(A)/(2|k0|d).  Our  results 
for  A/3  are  presented  in  Fig.  4c.  For  each  HHG  order,  we  estimate 
A/3  from  a  weighted  average  over  a  harmonic  peak,  without  any 
adjustable  fitting  parameters,  where  the  error  bars  are  the  standard 
deviation  over  the  same  harmonic.  The  MO  coefficient  measured  by 
the  left  circularly-polarized  3m  +  1  harmonics  matches  measure¬ 
ments  that  were  taken  using  synchrotrons48,  indicating  that  their 
polarization  is  circular  (see  Supplementary  Section  7  for  quantitat¬ 
ive  analysis  of  Stokes  parameter).  However,  the  low  intensity  and 
low  MO  dichroic  absorption  coefficient  measured  using  the  right 
circularly-polarized  29th,  32nd  and  35th  harmonics  (corresponding 
to  the  3m  -  1  group)  at  45  eV,  50  eV  and  55  eV,  respectively,  indi¬ 
cate  that  their  polarization  is  elliptical  with  an  ellipticity  of  -0.2. 
We  associate  the  asymmetry  between  the  right  and  left  high  harmo¬ 
nics  with  the  fact  that  the  rosette-shape  bi-chromatic  driving  laser 
induces  electronic  rotating  current  in  the  atoms,  influencing  both 
the  single-atom  HHG  process49  and  HHG  re-absorption. 
Interestingly,  the  XMCD  measurement  also  shows  that  the  polariz¬ 
ation  of  the  significantly  suppressed  3  m  harmonics  is  also  chiral, 
with  the  same  helicity  as  the  3m  +  1  harmonics.  This  observation 
agrees  with  the  numerical  calculations  presented  in  ref.  43 
(Figs  3c  and  3f),  where  the  polarizations  of  the  bi-chromatic 
drivers  deviate  from  perfect  circularity. 

Conclusions  and  outlook 

We  have  demonstrated  the  first  bright,  phase-matched  source  of  cir¬ 
cularly-polarized  high  harmonics  and  also  the  first  M-edge  X-ray 
magnetic  circular  dichroism  measurements  of  magnetic  materials 
on  a  table  top.  Our  HHG  source  produces  a  broad  spectrum  of 
bright  circularly-polarized  harmonics,  where  consecutive  harmo¬ 
nics  exhibit  opposite  helicity,  and  where  phase-matching  consider¬ 
ations  are  robust  and  favour  one  helicity  over  the  other.  This  work 
removes  a  major  constraint  to  date— that  the  polarization  of  bright 
high  harmonics  sources  is  limited  to  linear  polarization — and  there¬ 
fore  paves  the  way  for  ultrafast  circular  dichroism  studies  of 
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magnetism,  chiral  molecules  and  nanostructures,  as  well  as  for 
exploring  the  interaction  of  atoms  with  chiral  EUV  fields.  In  com¬ 
bination  with  coherent  diffraction  imaging  techniques,  this 
advancement  will  also  enable  the  exploration  of  coupled  spin, 
charge  and  structural  dynamics  of  magnetic  domains  with  elemen¬ 
tal  specificity  and  unprecedented  spatial  and  temporal  resolution. 
The  temporal  resolution  of  a  measurement  using  a  single  harmonic 
from  our  system  would  be  in  the  range  of  45  fs  (determined  by  the 
pump  and  the  harmonic  bandwidth).  To  access  the  attosecond 
regime  in  chiral  experiments,  interaction  with  several  helical  harmo¬ 
nics  that  rotate  in  the  same  direction  (and  therefore  form  a  helical 
attosecond  pulse  train)  is  required.  The  chirality- selective  phase - 
matching  can  very  nicely  enhance  only  harmonics  with  a  specific 
helicity,  thereby  producing  circularly-polarized  attosecond  pulses. 
Our  source  might  also  be  used  for  seeding  X-ray  free-electron 
lasers  to  obtain  ultra-bright  circularly-polarized  X-rays50.  Finally, 
the  scheme  we  use  is  universal,  and  can  be  used  to  generate  circu¬ 
larly-polarized  high  harmonics  across  very  broad  spectral  regions, 
for  example,  through  quasi  phase-matching  of  the  HHG  process, 
by  using  mid- infrared  driving  lasers15  or  by  using  other 
nonlinear  media. 

Received  18  April  2014;  accepted  30  October  2014; 
published  online  8  December  2014 

References 

1.  Bowering,  N.  et  al.  Asymmetry  in  photoelectron  emission  from  chiral 
molecules  induced  by  circularly  polarized  light.  Phys.  Rev.  Lett.  86, 

1187-1190  (2001). 

2.  Travnikova,  O.  et  al.  Circularly  polarized  X  rays:  another  probe  of  ultrafast 
molecular  decay  dynamics.  Phys.  Rev.  Lett.  105,  233001  (2010). 

3.  Gierz,  I.,  Lindroos,  M.,  Hochst,  H.,  Ast,  C.  R.  &  Kern,  K.  Graphene 
sublattice  symmetry  and  isospin  determined  by  circular  dichroism 
in  angle-resolved  photoemission  spectroscopy.  Nano  Lett.  12, 

3900-3904  (2012). 

4.  Xu,  S.-Y.  et  al.  Hedgehog  spin  texture  and  Berry’s  phase  tuning  in  a  magnetic 
topological  insulator.  Nature  Phys.  8,  616-622  (2012). 

5.  Liu,  Y.,  Bian,  G.,  Miller,  T.  &  Chiang,  T.-C.  Visualizing  electronic  chirality 
and  Berry  phases  in  graphene  systems  using  photoemission  with  circularly 
polarized  light.  Phys.  Rev.  Lett.  107,  166803  (2011). 

6.  Schiitz,  G.,  Kniille,  M.  &  Ebert,  H.  Magnetic  circular  X-ray  dichroism  and  its 
relation  to  local  moments.  Phys.  Scr.  1993,  302  (1993). 

7.  Stohr,  J.  et  al.  Element- specific  magnetic  microscopy  with  circularly  polarized 
X-rays.  Science  259,  658  (1993). 

8.  Fischer,  P.  et  al.  Magnetic  domain  imaging  with  a  transmission  X-ray 
microscope.  /.  Magn.  Magn.  Mater.  198-199,  624-627  (1999). 

9.  Eisebitt,  S.  et  al.  Lensless  imaging  of  magnetic  nanostructures  by  X-ray  spectro- 
holography.  Nature  432,  885-888  (2004). 

10.  Boeglin,  C.  et  al.  Distinguishing  the  ultrafast  dynamics  of  spin  and  orbital 
moments  in  solids.  Nature  465,  458-461  (2010). 

11.  Radu,  I.  et  al.  Transient  ferromagnetic-like  state  mediating  ultrafast  reversal  of 
antiferromagnetically  coupled  spins.  Nature  472,  205-208  (2011). 

12.  Beaurepaire,  E.,  Merle,  J.-C.,  Daunois,  A.  &  Bigot,  J.-Y.  Ultrafast  spin  dynamics 
in  ferromagnetic  nickel.  Phys.  Rev.  Lett.  76,  4250-4253  (1996). 

13.  Lopez-Flores,  V.  et  al.  Time-resolved  X-ray  magnetic  circular  dichroism  study  of 
ultrafast  demagnetization  in  a  CoPd  ferromagnetic  film  excited  by  circularly 
polarized  laser  pulse.  Phys.  Rev.  B  86,  014424  (2012). 

14.  Graves,  C.  E.  et  al.  Nanoscale  spin  reversal  by  non-local  angular  momentum 
transfer  following  ultrafast  laser  excitation  in  ferrimagnetic  GdFeCo. 

Nature  Mater.  12,  293-298  (2013). 

15.  Popmintchev,  T.  et  al.  Bright  coherent  ultrahigh  harmonics  in  the  keV 
X-ray  regime  from  mid-infrared  femtosecond  lasers.  Science  336, 

1287-1291  (2012). 

16.  Seaberg,  M.  D.  et  al.  Ultrahigh  22  nm  resolution  coherent  diffractive 
imaging  using  a  desktop  13  nm  high  harmonic  source.  Opt.  Express 
19,  22470-22479  (2011). 

17.  Kienberger,  R.  et  al.  Atomic  transient  recorder.  Nature  427,  817-821  (2004). 

18.  Worner,  H.  J.,  Bertrand,  J.  B.,  Kartashov,  D.  V.,  Corkum,  P.  B.  & 

Villeneuve,  D.  M.  Following  a  chemical  reaction  using  high-harmonic 
interferometry.  Nature  466,  604-607  (2010). 

19.  Siemens,  M.  E.  et  al.  Quasi-ballistic  thermal  transport  from  nanoscale 
interfaces  observed  using  ultrafast  coherent  soft  X-ray  beams.  Nature  Mater. 

9,  26-30  (2010). 

20.  Li,  W.  et  al.  Time-resolved  dynamics  in  N204  probed  using  high  harmonic 
generation.  Science  322,  1207-1211  (2008). 

104 


21.  Sandhu,  A.  S.  et  al.  Observing  the  creation  of  electronic  Feshbach 
resonances  in  soft  X-ray-induced  02  dissociation.  Science  322, 

1081-1085  (2008). 

22.  Cavalieri,  A.  L.  et  al.  Attosecond  spectroscopy  in  condensed  matter.  Nature 
449,  1029-1032  (2007). 

23.  Salieres,  P.  et  al.  Frequency-domain  interferometry  in  the  XUV  with  high-order 
harmonics.  Phys.  Rev.  Lett.  83,  5483-5486  (1999). 

24.  La-O-Vorakiat,  C.  et  al.  Ultrafast  demagnetization  dynamics  at  the  M  edges 
of  magnetic  elements  observed  using  a  tabletop  high-harmonic  soft  X-ray 
source.  Phys.  Rev.  Lett.  103,  257402  (2009). 

25.  Mathias,  S.  et  al.  Probing  the  timescale  of  the  exchange  interaction  in  a 
ferromagnetic  alloy.  Proc.  Natl  Acad.  Sci.  USA  109,  4792-4797  (2012). 

26.  Turgut,  E.  et  al.  Controlling  the  competition  between  optically  induced 
ultrafast  spin-flip  scattering  and  spin  transport  in  magnetic  multilayers.  Phys. 
Rev.  Lett.  110,  197201  (2013). 

27.  Corkum,  P.  B.  Plasma  perspective  on  strong  field  multiphoton  ionization. 

Phys.  Rev.  Lett.  71,  1994-1997  (1993). 

28.  Moller,  M.  et  al.  Dependence  of  high-order-harmonic-generation  yield  on 
driving-laser  ellipticity.  Phys.  Rev.  A  86,  011401  (2012). 

29.  Weihe,  F.  A.  et  al.  Polarization  of  high-intensity  high-harmonic  generation. 
Phys.  Rev.  A  51,  R3433-R3436  (1995). 

30.  Alon,  O.  E.,  Averbukh,  V.  &  Moiseyev,  N.  Selection  rules  for  the  high 
harmonic  generation  spectra.  Phys.  Rev.  Lett.  80,  3743-3746  (1998). 

31.  Yuan,  K.-J.  &  Bandrauk,  A.  D.  Circularly  polarized  molecular  high-order 
harmonic  generation  in  H2  with  intense  laser  pulses  and  static  fields.  Phys.  Rev. 
A  83,  063422  (2011). 

32.  Mauger,  F.,  Bandrauk,  A.  D.,  Kamor,  A.,  Uzer,  T.  &  Chandre,  C.  Quantum- 
classical  correspondence  in  circularly  polarized  high  harmonic  generation. 

/.  Phys.  B  Atom.  Mol.  Opt.  Phys.  47,  041001  (2014). 

33.  Vodungbo,  B.  et  al.  Polarization  control  of  high  order  harmonics  in  the  EUV 
photon  energy  range.  Opt.  Express  19,  4346-4356  (2011). 

34.  Liu,  L.  Z.,  O’Keeffe,  K.  &  Hooker,  S.  M.  Optical  rotation  quasi-phase-matching 
for  circularly  polarized  high  harmonic  generation.  Opt.  Lett.  37, 

2415-2417  (2012). 

35.  Fleischer,  A.,  Sidorenko,  P.  &  Cohen,  O.  Generation  of  high-order 
harmonics  with  controllable  elliptical  polarization.  Opt.  Lett.  38, 

223-225  (2013). 

36.  Yuan,  K.-J.  &  Bandrauk,  A.  D.  Single  circularly  polarized  attosecond  pulse 
generation  by  intense  few  cycle  elliptically  polarized  laser  pulses  and  terahertz 
fields  from  molecular  media.  Phys.  Rev.  Lett.  110,  023003  (2013). 

37.  Milosevic,  D.  B.,  Becker,  W.  &  Kopold,  R.  Generation  of  circularly  polarized 
high-order  harmonics  by  two-color  coplanar  field  mixing.  Phys.  Rev.  A 

61,  063403  (2000). 

38.  Milosevic,  D.  B.  &  Becker,  W.  Attosecond  pulse  trains  with  unusual  nonlinear 
polarization.  Phys.  Rev.  A  62,  011403  (2000). 

39.  Bandrauk,  A.  D.  &  Lu,  H.  Controlling  harmonic  generation  in  molecules 
with  intense  laser  and  static  magnetic  fields:  orientation  effects.  Phys.  Rev.  A 
68,  043408  (2003). 

40.  Zuo,  T.  &  Bandrauk,  A.  D.  High-order  harmonic  generation  in  intense  laser  and 
magnetic  fields.  J.  Nonlinear  Opt.  Phys.  Mater.  04,  533-546  (1995). 

41.  Long,  S.,  Becker,  W.  &  Mclver,  J.  K.  Model  calculations  of  polarization- 
dependent  two-color  high-harmonic  generation.  Phys.  Rev.  A  52, 

2262-2278  (1995). 

42.  Eichmann,  H.  et  al.  Polarization-dependent  high-order  two-color  mixing. 

Phys.  Rev.  A  51,  R3414-R3417  (1995). 

43.  Fleischer,  A.,  Kfir,  O.,  Diskin,  T.,  Sidorenko,  P.  &  Cohen,  O.  Spin  angular 
momentum  and  tunable  polarization  in  high-harmonic  generation.  Nature 
Photon.  8,  543-549  (2014). 

44.  Cohen,  O.,  Popmintchev,  T.,  Gaudiosi,  D.  M.,  Murnane,  M.  M.  &  Kapteyn,  H.  C. 
Unified  microscopic-macroscopic  formulation  of  high-order  difference- 
frequency  mixing  in  plasmas.  Phys.  Rev.  Lett.  98,  043903  (2007). 

45.  Popmintchev,  T.  et  al.  Phase  matching  of  high  harmonic  generation  in  the  soft 
and  hard  X-ray  regions  of  the  spectrum.  Proc.  Natl  Acad.  Sci.  USA  106, 
10516-10521  (2009). 

46.  Oppeneer,  P.  M.  in  Handbook  of  Magnetic  Materials  Vol.  13  (ed.  Buschow,  K.  H.  J.) 
229-422  (Elsevier,  2001). 

47.  Henke,  B.  L.,  Gullikson,  E.  M.  &  Davis,  J.  C.  X-ray  interactions:  photoabsorption, 
scattering,  transmission,  and  reflection  at  E  -  50-30,000  eV,  Z  =  1-92.  Atom. 
Data  Nucl.  Data  Tables  54,  181-342  (1993). 

48.  Valencia,  S.  et  al.  Faraday  rotation  spectra  at  shallow  core  levels:  3 p  edges  of  Fe, 
Co,  and  Ni.  New  J.  Phys.  8,  254  (2006). 

49.  Xie,  X.  et  al.  Internal  momentum  state  mapping  using  high  harmonic  radiation. 
Phys.  Rev.  Lett.  101,  033901  (2008). 

50.  Zeitoun,  P.  et  al.  High-intensity  highly  coherent  soft  X-ray  femtosecond  laser 
seeded  by  a  high  harmonic  beam.  Nature  431,  426-429  (2004). 

Acknowledgements 

This  work  was  supported  by  the  USA-Israel  Binational  Science  Foundation  (BSF).  The 

Technion  group  was  supported  by  the  Israel  Science  Foundation  (grant  no.  1225/14)  and 

NATURE  PHOTONICS  |  VOL  9  |  FEBRUARY  2015  |  www.nature.com/naturephotonics 


©  2015  Macmillan  Publishers  Limited.  All  rights  reserved 


NATURE  PHOTONICS  DOI:  10.1038/NPHOTON.2014.293 


ARTICLES 


the  Israeli  Center  of  Research  Excellence  ‘circle  of  light’  supported  by  the  I-CORE  Program 
of  the  Planning  and  Budgeting  Committee  and  The  Israel  Science  Foundation.  The 
experiment  was  carried  out  in  JILA.  O.K.  and  O.C  acknowledge  the  warm  hospitality  in 
JILA.  The  JILA  and  NIST  authors  acknowledge  funding  from  the  US  Department  of  Energy 
Office  of  Basic  Energy  Sciences  (award  no.  DE-SC0002002).  JILA  also  acknowledges 
funding  from  the  Physics  Frontiers  Center  Program  and  from  an  AFOSR  DURIP  award  for 
the  laser  system  used  for  this  work.  P.G.  acknowledges  support  from  the  Deutsche 
Forschungsgemeinschaft  (no.  GR  4234/1-1).  R.K.  acknowledges  the  Swedish  Research 
Council  (VR)  for  financial  support. 

Author  contributions 

O.C.,  M.M.  and  H.K.  conceived  and  supervised  the  project.  O.K.,  P.G.,  E.T.  and  A.F. 
designed  the  experimental  set-up.  O.K.,  P.G.,  E.T.,  R.K.,  D.Z.,  T.P.  and  D.P.  constructed  the 


experimental  set-up  and  conducted  the  experiment.  O.K.  and  P.G.  analysed  the  data  with 
input  from  other  authors.  O.K.  developed  the  theory  and  carried  out  the  numerical 
simulations.  O.K.,  P.G.  and  O.C.  wrote  the  manuscript  with  input  from  all  other  authors. 

Additional  information 

Supplementary  information  is  available  in  the  online  version  of  the  paper.  Reprints  and 
permissions  information  is  available  online  at  www.nature.com/reprints.  Correspondence  and 
requests  for  materials  should  be  addressed  to  O.K.  and  O.C. 

Competing  financial  interests 

The  authors  declare  no  competing  financial  interests. 


NATURE  PHOTONICS  |  VOL  9  |  FEBRUARY  2015  |  www.nature.com/naturephotonics 


105 


©  2015  Macmillan  Publishers  Limited.  All  rights  reserved 


photonics 


nature 


Generation  of  bright  phase-matched 
circularly-polarized  extreme  ultraviolet 

high  harmonics 


Ofer  Kfir1,  Patrik  Grychtol2,  Emrah  Turgut2,  Ronny  Knut2,3,  Dmitriy  Zusin2, 

Dimitar  Popmintchev2,  Tenio  Popmintchev2,  Hans  Nembach2,3,  Justin  M.  Shaw3,  Avner  Fleischer1,4, 
Henry  Kapteyn2,  Margaret  Mumane2  and  Oren  Cohen1. 

1 Solid  State  Institute  and  Physics  Department,  Technion,  Haifa  32000,  Israel 
2 Department  of  Physics  and  JILA,  University  of  Colorado  and  NIST,  Boulder,  CO  80309,  USA 
3  Electromagnetics  Division,  National  Institute  of  Standards  and  Technology,  Boulder,  CO  80305,  USA 
4 Department  of  Physics  and  Optical  Engineering,  Ort  Braude  College,  Karmiel  21982,  Israel 

1.  Phase  matching  of  circular  HHG 

Full  phase  matching  of  a  nonlinear  process  corresponds  to  where  the  coherence  length  of  the 
process  diverges;  i.e.  when  the  right  hand  side  of  Eq.  (3)  equals  zero.  To  make  this  possible, 
the  effective  index  changes,  An1  and  A n2,  must  have  opposite  sign.  Furthermore,  only  one 
harmonic  order  can  be  fully  phase  matched  for  a  given  effective  index  ratio  of  the  two 
colors  An2/An1.  In  order  to  phase  match  the  harmonic  order  q  =  3m  +  1,  the  effective  index 
ratio  should  be  - 


A n2  _  m±  1 
A  n1  2m 


(S.l) 


For  example,  in  order  to  phase  match  the  28th  or  40th  HHG  orders,  (m  =  9,13)  and  (+— >  +), 
the  index  ratio  An2/An1  should  be  -0.55  or  -0.54,  respectively. 

Explicit  representation  of  the  effective  index  change  for  each  component  of  the  bi-chromatic 
driver,  An^  A n2  is  given  by1 


Here,  i=  1  for  the  fundamental  and  1=2  for  the  second  harmonic  driver,  u11  is  the  modal  factor2, 


a  is  the  inner  radius  of  the  hollow  waveguide,  and  kt  are  the  wavelength  and  wave-vector 
of  z'th  field,  P  is  the  gas  pressure  in  atmospheres,  17  is  the  ionization  fraction,  re  is  the  classical 
radius  of  the  electron,  Natm  is  the  number  density  of  atoms  at  atmospheric  pressure,  and  rq  is 
the  index  of  refraction  of  the  gas  for  wavelength  /  at  atmospheric  pressure.  The  modal  term, 
ufjAf /Sn2a2,  and  the  plasma  term  r}Natmre  / An,  contribute  negatively  to  An,  mainly  to  the 
fundamental  driver,  Al9  while  the  normal  refraction  of  the  neutral  gas,  (n,  —  1),  contributes 
positively  to  An,  mainly  to  the  second  harmonic  field,  A2.  Since  full  phase  matching  requires 
that  An2/An1  <0,  it  is  achievable  when  Anx  (long  wavelength)  is  negative  and  A n2  (short 
wavelength)  is  positive. 
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The  gas  pressure  in  which  high  harmonics  are  generated  acts  as  a  continuous  knob  to  control 
the  HHG  phase  matching.  At  zero  pressure,  the  effective  index  ratio  is  A n2  /  A%  =  1/4,  where 
both  An2,An1  are  negative  due  to  modal  term,  u^Xf /8n2a2  .  Pressurizing  the  hollow 
waveguide  with  a  weakly  ionized  gas  contributes  positively  to  A n2,  more  than  to  An1;  so  high 
enough  gas  pressure  can  make  An2/An1  <  —0.5,  as  required  for  phase  matching  of  high  order 
harmonics.  Tuning  the  pressure,  therefore,  selects  the  HHG  order  to  be  phase  matched 
according  to  Eq.  (S.l)  and  (S.2).  However,  re-absorption  of  HHG  in  the  gas  reduces  the  contrast 
between  the  optimally  phase  matched  harmonic  (e.g.  28th  order  in  Fig.  lb)  and  other  harmonic 
orders. 

Figure  S.l  plots  the  coherence  length  vs.  the  absorption  length  of  few  HHG  orders  in  Ne  for  a 
1%  ionization  level.  When  the  ratio  of  coherence  length,  Lc,  to  absorption  length,  La,  is  above 
6  (see  dashed  line),  the  HHG  flux  is  90%  of  the  fully  phase  matched  case,  and  the  emission  is 
effectively  phase  matched.  The  low  order  HHG  (19th  and  20th)  increase  at  pressures  around  100 
Torr  and  are  phase  matched  up  to  far  above  1000  Torr,  whereas  the  high  order  harmonics  (40th 
and  41st)  are  phase  matched  in  a  narrower  pressure  range.  The  pressure  dependent  spectrograms 
in  Fig.  3c  (measured)  and  Fig.  3d  (theory)  agree  with  this  description,  although  the  calculation 
takes  into  account  additional  effects,  such  as  the  HHG  source  density  (i.e.  gas  density)  and 
parasitic  re-absorption  in  residual  gas  downstream  from  the  phase  matching  region  (details 
below). 


Pressure  [Torr] 

Figure  S.l  I  (a)  The  coherence  length  vs.  absorption  length  of  circularly  polarized  harmonics 
in  a  Ne  filled  hollow  waveguide.  At  Lc/La  =  6  (see  dashed  black  line),  the  output  HHG  flux 
is  90%  of  the  flux  under  perfect  phase  matching  condition,  so  the  HHG  are  said  to  be  phase 
matched.  The  low  order  HHG  orders  (say  19th  and  20th  )  are  phase  matched  from  about  100 
Torr  to  above  1000  Torr,  whereas  the  high  order  harmonics  (40th  and  41st)  are  phase  matched 
in  a  narrower  pressure  range.  The  solid  and  fine  lines  correspond  to  HHG  of  the  groups  q  = 
3m  +  1  and  q  =  3m  —  1,  respectively.  Lc  of  the  19th  order  deviates  from  the  20th  due  to  the 
helicity  selective  factor,  ±An,  in  Eq.  (3). 
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2.  Critical  ionization  for  full  phase  matching  of  circularly  polarized  HHG 

The  critical  ionization,  rjcr,  is  the  ionization  fraction  above  which  pressure  tuning  cannot  bring 
about  phase  matching  of  circularly  polarized  HHG.  Specifically,  pressure  tuning  cannot  reduce 
the  ratio  An2/An1  below  -0.5.  The  value  of  t]cr  can  be  derived  from  Eq.  (S.2)  - 

T]cr  =  [2  (n2  -  1)  +  (%  -  1)]  •  {2(n2  -  1)  +  (%  -  1)  +  ^  Natmre  Af}  (S.3) 

This  critical  ionization  differs  from  critical  ionization  of  HHG  driven  with  single  laser  field2. 
For  HHG  driven  by  one  field  in  a  large  focus  geometry,  the  critical  ionization  level  is 
independent  of  the  pressure,  and  occurs  when  the  dispersion  of  the  neutral  atoms  and  plasma 
balance  each  other.  (In  a  waveguide,  the  value  of  the  critical  ionization  is  reduced  slightly 
below  that  of  a  large  focal  geometry,  due  to  the  presence  of  waveguide  dispersion).  In  contrast, 
circularly  polarized  HHG  driven  by  bi-chromatic  fields  requires  that  An*  is  nonzero:  A% 
should  be  negative  and  A n2  positive.  Since  optimal  A^  is  negative  and  A n2  is  positive,  r]cr  is 
larger  for  bi-chromatic  circular  driver,  than  for  HHG  driven  with  the  fundamental  only,  but 
smaller  then  HHG  driven  by  the  second  harmonic  only. 

3.  Numerical  procedure  for  a  pressure  dependent  spectrogram 

Simulating  the  single  atom  response  to  the  driving  bi-chromatic  field  is  the  principal  step  for 
calculating  pressure  dependent  HHG  emission.  We  approximate  the  atom  as  a  binding  potential 
and  a  single,  spin-less  electron,  and  solve  the  Schrodinger  equation  in  three  spatial  dimensions 
(3D)  and  time.  The  binding  potential  is  an  inverted  Gaussian  given  by  Vatom(r )  =  —  V0e~s^r^2 , 
where  V0  =  2.38  and  s  =  1/V8  atomic  units  (at.u.).  The  single  electron  initially  occupies  the 
lowest  state,  having  ionization  potential  of  a  Ne  atom,  0.794  at.u.  (21.6  eV).  The  driving  bi- 
chromatic  field,  EBC(t),  induces  a  linear  potential,  within  the  dipole  approximation,  so  overall 
the  potential  is  V(r,t)  =  - EBC(t )  •  r  +  Vatom(r)  .  Here  EBC(t)  =  Et(t)  +  E2(t)  ,Et(t)  = 

(— sin(a)0t)  x  +  cos(co0t)y)  and  E2(t)  =  E2, oe-^?)  (sin(2o)0t  +  02)  x  + 
cos(2 co0t  +  4>2)y),  where  E10,  E2  0  correspond  to  peak  intensities  of  2.3  •  10 18 [W /m2]  and 
2.3  •  10 18[W/m2],  respectively.  tp,  correspond  to  15  femtoseconds,  o>0  =  0.0577  is  the 
frequency  of  790nm  light  in  atomic  units,  02  is  the  phase  of  the  second  harmonic  field  relative 
to  the  fundamental,  x,  y  are  the  Cartesian  coordinates  unit  vectors  and  f  is  the  radius  vector. 
The  6-power  Gaussian  envelope  of  the  simulated  bi-chromatic  pulse  is  chosen  to  reduce 
spectral  effects  originating  from  variations  of  the  pulse  envelope.  The  emitted  HHG  field, 
Ehhg,  is  calculated  using  Ehrenfest  theorem,  EHHG(t)  oc  (0(t)|  Vl/ato?n|0(t)},  where  0  is  the 
electronic  wave  function.  The  proportion  constant  is  not  accounted  for  since  it  applies  globally 
to  the  entire  radiation  spectrum. 

The  HHG  phase  strongly  depends  on  the  effective  index  of  the  bi-chromatic  driving  field  in 
the  gas  medium.  To  account  for  the  phase  differences  between  the  fundamental  and  second 
harmonic  fields,  we  calculate  HHG  emission  from  a  single  atom  100  times,  where  each 
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iteration  uses  a  different  phase  between  the  fundamental  and  second  harmonic  drivers.  The 
phase  for  the  hth  iteration  is  02,h  =  2nh/100,  where  h=l,2,3...100,  and  the  calculated  HHG 
emission  is  assigned  to  an  atoms  positioned  at  zh  =  <p2ihC[2co0(An1  —  An2)]_1.  The  effect  of 
the  fundamental  driver  phase  velocity  is  accounted  for  once  the  single  atom 
emission,  EHHGh(t),  is  done.  Since  the  fundamental  phase  velocity  in  the  gas  filled  waveguide 
(Eq.  (S.2))  differs  from  free  space  propagation,  it  requires  additional  time,  An1zh/c ,  to  reach 
the  atom  at  zh.  Therefore,  the  emission  at  zh  is  delayed  from  the  calculated  emission,  and  is 
given  by  EHHG,h(t  -  An^Jc). 

The  macroscopic  HHG  phase  matched  signal  is  a  coherent  sum  of  the  single  atom  emissions 
that  undergo  partial  re-absorption  in  the  gas.  At  the  end  of  the  phase  matching  region,  zPM,  the 

HHG  intensity  is  [ZhPEHHGh(t  —  An1zh/c)e~(zpM~Zh^Lahs\  ,  where  Labs  is  the  absorption 
length  of  the  HHG  in  the  gas.  The  gas  pressure,  P,  is  inserted  to  the  calculation  to  account  for 
the  density  of  HHG  emitters.  Upon  exiting  the  phase  matching  region,  the  HHG  field  can 
undergo  additional  parasitic  re-absorption  in  the  neutral  gas.  This  is  simulated  by  adding 
absorption  from  a  150  pm  layer  of  the  gas  at  pressure  P.  The  pressure-dependent  spectrogram 
in  Fig.  3d  is  calculated  by  repeating  the  above  for  multiple  pressure  of  Ne  gas  at  an  ionization 
fraction  rj  =  1.2%. 


4.  Dynamics  of  circularly  polarized,  counter  rotating  bi-chromatic  field 

The  polarization  of  the  light  field  is  usually  described  by  its  Lissajous  curve,  that  is,  the  electric 
field  of  circular,  elliptical,  and  linear  polarization  follows  a  circle,  ellipse  and  a  line, 
respectively.  In  that  sense,  the  polarization  of  a  circularly  polarized,  counter  rotating  bi- 
chromatic  field  composed  of  a  fundamental  field  at  790  nm  and  its  second  harmonic  has  a  3- 
fold  “rosette”  shape  polarization,  as  described  in  Fig.  la  in  the  text.  To  better  illustrate  the 
dynamical  aspect  of  this  rosette  polarization  we  attach  a  movie  showing  the  fundamental  driver 
(red)  second  harmonic  driver  (blue)  and  the  resulting  field  of  the  bi-chromatic  driver  (purple). 
(See  supplementary  movie ) 

5.  Calculating  and  comparing  the  photon  flux  of  circular  and  linear  polarized  harmonics. 

To  determine  the  feasibility  of  future  XMCD  pump-probe  experiments  using  a  source  of 
circularly  polarized  high  harmonics,  we  have  calculated  the  absolute  flux  contained  in  a  single 
circular-polarized  harmonic  at  the  Co  Af-edge  below  60  eV  and  compared  it  with  the  flux 
obtained  in  previous  ultrafast  element-selective  HHG  experiments  exploiting  the  transversal 
magneto-optical  Kerr  effect  (T-MOKE)  in  reflection  using  linear  polarized  harmonics3  7.  To 
this  end,  we  have  simulated  the  transmission  as  well  as  reflection  properties  of  our  EUV 
beamline  components  taking  into  account  the  efficiency  of  our  CCD  camera.  Specifically,  we 
have  used  the  RAY  simulation  software  package  that  is  a  powerful  design  tool  developed  at 
BESSY  for  synchrotron  radiation  beamlines8.  Moreover,  we  have  calibrated  the  efficiency  of 
our  CCD  detector  according  to  the  procedure  recommended  by  the  manufacturer. 
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As  schematically  depicted  in  Fig.  2  of  the  main  article,  the  harmonics  emerging  from  our 
waveguide  pass  through  a  200  nm  thick  A1  filter  (with  a  5  nm  thick  AI2O3  layer)  and  the  100 
nm  thick  Co  film,  which  is  tilted  by  45°  with  respect  to  the  HHG  beam,  before  entering  our 
spectrometer.  The  spectrometer  consists  of  a  toroid  mirror  (glass  substrate  coated  by  100  nm 
of  B4C  and  tilted  at  a  grazing  angle  of  8°),  a  laminar  gold  grating  patterned  on  top  of  Si  substrate 
(with  a  groove  density  of  500  lines/mm,  a  groove  depth  of  20  nm  and  which  is  mounted  at  a 
grazing  angle  of  12°),  two  200  nm  thick  A1  filter  (each  of  which  having  a  5  nm  thick  AI2O3 
layer)  and  a  CCD  camera  (Andor  -  Newton  DO-920).  As  can  been  seen  in  Fig.  S.2,  this 
combination  of  components  results  in  a  theoretical  throughput  of  10'5-10"4  EUV  photons  in  the 
vicinity  of  the  Co  M-edge  at  58  eV.  This  is  comparable  to  the  efficiency  of  previously  reported 
measurements  in  the  T-MOKE  geometry. 


photon  energy  (eV) 

Figure  S.2  I  Efficiency  comparison  of  T-MOKE  and  XMCD  HHG  EUV  beamlines 


In  this  spectral  region,  the  CCD  chip  of  the  camera  has  a  quantum  efficiency  of  about  30%,  i.e. 
about  every  third  photon  is  absorbed  by  the  back-thinned  Si  chip.  Each  absorbed  photon  with 
an  energy  of  55eV  generates  approximately  55/3.65  photoelectrons  or  55/45.63  counts  using 
the  standard  gain  setting  of  our  camera.  Taking  into  account  the  simulated  efficiency  of  our 
XMCD  (T-MOKE)  beamline  of  about  4T0"5  (1.5*10~5)  and  the  fact  that  we  observe  around  104 
counts/s  in  the  55  eV  harmonics  of  both  HHG  setups,  we  obtain  a  photon  flux  of  approximately 
7108  photons/s  in  a  single  circular  harmonic,  which  is  only  2-3  times  smaller  than  the  flux  of 
the  very  same  linear  harmonic  in  our  T-MOKE  experiment.  Considering  the  very  high  signal  - 
to-noise  ratio  that  we  have  already  attained  in  our  static  XMCD  experiment  as  well  as  the 
relatively  high  effective  thickness  of  our  Co  sample  (140  nm)  in  comparison  to  the  penetration 
depth  of  a  laser  excitation  pulse  (a  few  tens  of  nm),  we  conclude  that  ultrafast  pump-probe 
investigations  of  chiral  phenomena  are  within  the  capabilities  of  our  current  experimental 
circular  HHG  apparatus. 
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6.  Useful  features  of  circular  HHG  using  a  counter-rotating  bi-chromatic  field. 

We  generate  bright  circularly  polarized  harmonics  by  co-propagating  bi-chromatic  790  nm  and 
395  nm  (to  and  2to)  driving  laser  beams  with  counter-rotating  circular  polarization9-11.  This 
approach  has  several  useful  properties  that  are  beneficial  for  applications.  First,  the  conversion 
efficiency  to  circularly  polarized  high  harmonics  is  comparable  to  the  conversion  efficiency  of 
traditional  HHG,  where  linearly  polarized  driving  lasers  produce  linearly  polarized  high 
harmonics  (see  Fig.  S.2  in  Ref11).  Second,  the  generated  spectrum  consists  only  of  circularly 
polarized  harmonics.  Third,  consecutive  harmonics  exhibit  opposite  helicity.  Fourth,  the 
circular-polarization  of  the  harmonics  is  manifested  in  the  high  harmonic  spectrum  due  to 
angular  momentum  conservation  rules:  thus,  it  is  not  essential  to  directly  or  continuously 
monitor  their  polarization.  As  explained  below,  when  the  bi-chromatic  driving  lasers  consist  of 
the  fundamental  and  second  harmonic  frequencies,  absence  of  the  3mth  harmonics  (m=l  ,2,3 . . .) 
indicates  that  the  polarizations  of  3m+l  and  3m- 1  harmonics  are  circular  with  helicity  that 
correspond  to  the  polarization  of  the  fundamental  and  2nd  harmonic  fields,  respectively.  Fifth, 
the  polarization  of  the  generated  harmonics  is  largely  insensitive  to  the  intensities  and  intensity 
ratio  between  the  bi-chromatic  driving  lasers.  This  is  an  important  feature  because  both  the 
intensity  of,  and  the  intensity  ratio  between,  the  two  driving  fields  can  vary  within  the  focal 
spot  and  along  the  direction  of  propagation.  Also,  this  feature  is  important  for  the  robustness 
of  the  method  to  misalignment  and  instability  of  the  driving  fields.  Sixth,  as  shown  below,  the 
generation  process  can  be  fully  phase  matched.  Finally,  the  mechanism  for  producing  circularly 
polarized  harmonics  in  this  direct  approach  is  based  only  on  the  circular  polarization  of  the 
driving  pulses  and  is  otherwise  insensitive  with  respect  to  the  nonlinear  medium  and  other 
properties  of  the  driving  laser  beams.  Thus,  we  expect  that  many  experimental  applications  that 
were  implemented  using  linearly  polarized  HHG  can  be  directly  transferred  to  circularly 
polarized  harmonics.  For  example,  we  expect  that  the  use  of  long-wavelength  bi-chromatic 
driving  lasers  will  extend  bright  circularly-polarized  HHG  into  the  keV  region,  which  will  span 
the  water-window  region  as  well  as  the  magnetically  sensitive  L-shell  absorption  edges  of  the 
3d  ferromagnets12. 


7.  Degree  of  circular  polarization 

The  4  Stocks  parameters13  S0,  S±,  S2,  S3  fully  characterize  the  polarization  state  of  light: 
the  0th  parameter  S0  measures  the  overall  radiation  intensity,  S1  and  S2  measure  the  degree  of 
linear  polarization  and  its  orientation  and  S3,  measures  the  degree  of  circular  polarization  and 
helicity.  More  accurately,  the  degree  of  circular  polarization  is  given  by  S3  normalized  by  the 
total  intensity  S0  and  has  the  form 


So 


Ir—Il 

lR+lL 


(S.4) 


where  lR  is  the  intensity  of  the  right-rotating  circularly-polarized  component  of  the  light  and  IL 
is  the  intensity  of  the  left-rotating  circularly-polarized  component.  Usually  one  records  IR,  IL 
by  measuring  the  intensity  passing  through  a  polarizer  that  transmits  pure  circular  polarization 
with  right  or  left  helicity,  respectively.  Measuring  S3/S0  =  1  indicates  pure  right-rotating 
circularly-polarized  light  (RCP)  while  measuring  S3/S0  =  —1  corresponds  to  pure  left- 
rotating  circularly-polarized  light  (LCP). 
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The  XMCD  measurement  that  we  conducted  within  this  work  retrieves  the  degree  of  circular 
polarization  in  the  spectral  region  where  Co  is  active  magneto-optically.  Incident  light  on  the 
Co  foil  can  be  decomposed  into  RCP  light  of  intensity  IR  and  LCP  light  of  intensity  IL.  The 
light  intensity  transmitted  through  the  Co  sample  with  “up”  or  “down”  magnetization  depends 
on  the  MO  constant  A/?14,  the  light  wave  vector  k0  and  the  sample  thickness  L 

J  up  —  jRe2\k0\A/SL  _|_  ^g-2|fc0|A/?L 
j  down  _  jRg-2|fc0|A/?r  _|_  ^g+2|fc0|A/?L 

Substituting  eq.  S.5  into  the  definition  of  the  XMCD  asymmetry, 

/down)  ?  an(j  using  eq.  S.4,  gives  the  degree  of  circular  polarization: 

S3  _  A 

S0  tanh(2|fc0|A/?L)' 

Fig  S.3  presents  the  degree  of  circularity  retrieved  from  the  measured  asymmetry.  The  left 
rotating  HHG  orders  31st  (48  eV)  and  34th  (53  eV)  are  highly  circular  with  S3/S0  =  —1.01  ± 
0.15  and  -0.98  ±  0.19  ,  respectively.  That  is,  within  our  experimental  uncertainty,  the 
polarization  is  pure  LCP!  On  the  other  hand,  the  right  rotating  HHG  orders  32nd  (50  eV)  and 
34th  (55  eV)  have  low  degree  of  circularity,  S3/S0  —  0.37  +  0.22  and  0.31  ±  0.10,  respectively. 
The  degree  of  circularity,  S3/S0,  cannot  be  retrieved  below  45  eV  since  A/?  is  still  unknown14. 
Furthermore,  at  the  37th  HHG  order  (58  eV)  the  MO  coefficient,  A/?,  nullifies  so  eq.  S.6 
diverges  and  S3/S0  cannot  be  retrieved. 


(S.5) 


=  (/up  _  /down) /(/up  + 


(S.6) 
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Figure  S.3  I  Normalized  Stocks  parameter  S3/S0  for  HHG  in  the  Cobalt  M-edge.  The 

magneto  optical  selectiveness  in  the  M-edge  of  Co  acts  essentially  as  a  polarizer  for  circular 
polarization.  Using  the  XMCD  asymmetry  we  conclude  that  HHG  orders  q  =  3 m  +  1  are 
purely  left-rotating  circularly  polarized  (within  our  experimental  uncertainty):  S3/S0  —  —1.01  + 
0.15  for  the  31st  order  (48  eV)  and  S3/S0  =  -0.98  +  0.19  for  the  34th  order  (53  eV).  The  right 
rotating  HHG  orders  orders,  q  =  3 m  —  1,  present  low  degree  of  circularity:  S3/S0  =  0.37  +  0.22 
for  the  32nd  order  (50  eV)  and  S3/S0  =  0.31  ±  0.10  for  the  35th  order  (55  eV).  Therefore,  the 
right  rotating  HHG  are  elliptical  or  partially  polarized.  The  degree  of  circularity  and  the  error 
bars  are  calculated  from  the  average  and  standard  deviation  of  the  XMCD  asymmetry  over 
each  harmonic  order  divided  by  tanh(2|/e0|A/?L)  according  to  eq.  S.6. 


Exp.  -  q  =  3m  +  1 
Exp.  -  q  —  3 m  -  1 


O 


O 


+  +■  , 

Right  rotating  elliptical  HHG 

Left  rotating 

circularly  polarized  HHG 
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